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ABSTRACT 



The transient testing technique for the determination of heat trans- 
fer characteristics of a porous medium has been known since 1929. It 
requires a good deal of care and one of the important restrictions is a 
step change in temperature must be applied to the medium at the beginning 
of a test run. This thesis involves the design, test and evaluation of 
an experimental set-up that utilizes an "integral heater" to accomplish 
the step change in temperature. In addition, there is presented heat 
transfer and flow friction characteristics for a 20° skew matrix con- 
structed from perforated nickel material. This data is compared to a 
20° skew matrix made from stainless steel material that has been pre- 
viously reported by C. P. Howard. Also, data is presented for two "paral- 
lel plate" matrices, one constructed from brass sheet, the other from 
perforated nickel material. A comparison of their relative performances 
is made. 
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Matrix Flow length, ft 
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pressure, Ib/hr 

Hydraulic radius, ft 
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Matrix Volume ft 
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w 

s 


= 


Matrix mass, lb 


(3 


= 


2 3 

Matrix compactness, ft /ft 


(3 




Compactness corrected for effect of slots 


A- 


- 


Fluid viscosity, lb/hr ft 




= 


Heat transfer area reduction factor 


s 


= 


Conduction area reduction factor 


f 


- 


Fluid density, lb/ft ^ 




- 


3 

Matrix density, lb/ft 


Nondimens Iona 1 parameters: 


f 


- 


Fanning friction factor, ratio of the wall shear stress 
to the fluid dynamic head 


j 




2/3 

Colburn j factor, N g T N p^ 


n nu 


= 


Nusselt's number, h D^, 




k 


n r 


= 


Reynold's number 4r H G/ 


n st 


= 


Stanton number, h/Gc 

P 


n pr 


= 


Prandtl number, IX c /k 
P 


N tu 


= 


Number of heat transfer units, hA/mc 

P 


K 

e 


= 


Entrance flow coefficient for the matrix 


K c 


= 


Exit flow coefficient for the matrix 


P 


= 


porosity of the matrix 

Matrix flow void volume/Matrix volume 


*X 


- 


Conduction parameter, k g A g /LmCp 


Subscript 


s : 


i 


= 


initial 


1 


= 


upstream 


2 


= 


downstream 


f 


* 


fluid 
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STD = 



solid (matrix material) 

Standard temperature and pressure 
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INTRODUCTION. 



1 . 

A compact heat exchanger is one not necessarily of small bulk and 
weight, but one which incorporates a heat transfer surface having a high 
area density. 

Arbitrarily, it has been specified that $ — 200 ft^/ft"* makes 

for a compact heat exchanger, where >0 is the ratio of the heat trans- 
2 3 + 

fer area to core volume, ft /ft . (15) 

Motivation for using compact surfaces is to gain high exchanger 
effectiveness in given box volume and weight limitations and at a reason- 
able cost. Great impetus has been given the compact heat exchanger field 
by recent applications as a regenerator in aircraft turboprop power plants, 
as a regenerator in the gas turbine automobile engine and by their use 
in cryogenic plant exchangers in the production of liquid oxygen and 
nitrogen. As the applications for compact heat exchangers increase, 
design problems increase, A successful compact heat exchanger design 
and development program requires access to, or the ability to generate 
knowledge of the basic heat transfer and friction behavior of compact 
surfaces of interest. It is also desirable that the designer would have 
available as much information on surface performances as possible to re- 
main flexible in his design with regard to the limitations that have been 
placed upon him. 

Building a complete heat exchanger to evaluate a new surface as to 
its heat transfer and friction performance would be prohibitively ex- 
pensive. The alternate method is to generate the performance character- 
istics of a surface of interest by testing small matrix modules and use 

+Numbers in parentheses refer to references in the bibliography. 
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this data to design a prototype heat exchanger. 

Briefly, the objectives of this thesis are: to build an experimental 
rig that will test small matrix modules to determine their heat transfer 
and flow friction characteristics, and, to test matrices made from per- 
forated nickel material to determine these characteristics. 

One of the simplest and most satisfactory methods of measuring the 
forced convection heat transfer characteristics of a matrix utilizes a 
transient test technique or so called single blow method. This test 
technique consists of heating the homogeneous porous solid of known 
weight and geometry with a uniformly distributed heat source and then 
at a given instant, removing the heat source and simultaneously cooling 
the solid with a fluid at a uniform lower temperature. The heating and 
cooling steps can be reversed and the same results achieved. A time- 
temperature response curve of the fluid stream leaving the solid is 
recorded. By comparing this record to theoretically developed non-dimen- 
sional temperature response curves, the heat transfer characteristics of 
the solid can be determined. 



2 



2 . 



SUMMARY OF THEORY. 



A theoretical analysis of the transient heat transfer behavior of 
a porous medium was presented by Schumann (2) in 1929. A porous medium 
may be defined as an aggregation of solid bodies so arranged that the 
void spaces between the bodies are connected and will permit the passage 
of a fluid. This definition can of course be extended to include the 
matrix surfaces utilized in a compact heat exchanger. 

In Schumann's analysis, the porous medium is considered to be homo- 
geneous and initially at a uniform temperature, t^. A fluid at the same 
temperature is flowing through the medium. At a certain instant the 
temperature of the entering fluid is assumed to change to a higher tem- 
perature t f ^. The problem is then to find the temperature of the fluid 
and the solid as functions of time and position in the matrix or porous 
media. The basis of the analysis are the theoretical generalized heat- 
ing curves that are developed from an energy balance on the matrix. 





I 




Solid 
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Assume : 



1) Properties of the fluid are independent of temperature. 

2) Fluid flow is steady. 

3) The porous solid is homogeneous. 

4) The thermal conductivity of both fluid and solid is infinite 
perpendicular to the flow direction. 

5) Thermal conductivity of fluid is zero in the flow direction. 



Heat absorbed by the solid = ^ 4s Cs dQ 

Heat transferred to solid by convection = s)Jk 

. ?tf / 

Heat transferred from the fluid by convection = mCp wX 
Heat transferred from solid by conduction = " k s As 



An energy balance for the fluid and solid results ini 






mc p d* - t s ) Jk =0 



= L L 

z = generalized position variable = N 



Let: T= generalized time variable = pP © 

W S C S 



( 1 ) 

( 2 ) 



N 



tu 



h A 
m c 



X- 



P 



Making the above substitutions, equations 1) and 2) become: 

itf _ (±s-U) 

The experimental method used requires that the above assumptions 



and the following boundary and initial conditions be closely observed. 
These conditions are as follows: 

1) Initially the core is at a uniform temperature. 

2) At time t = 0, the temperature of the incoming fluid changes 
instantaneously. 
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3) System is adiabatic, that is, only heat is transferred be- 
tween the matrix core and the fluid, none with the surroundings. 

With the idealizations, initial and boundary conditions, and the additional 
assumption that the thermal conduction in the solid is zero parallel to 
the direction of fluid flow (i.e., X , the material conduction para- 
meter, = 0), Schumann's solutions are as follows: 



-(z +T ) 



CO 

2 ' 



d (z T ) 



(2i >Jt z) 




1 - 






+ T) 



CO 




n*o 



n 



z 



d n 

d(zT ) n 



j 0 (2ipr7) 



This solution was first used as the basis for a transient technique by 
Furnas (16) in 1932. The fluid temperature was measured at the exit of 
the matrix where x = L, therefore z = N^. The experimental data was 
compared with Schumann's theoretical constant z curves. The theoretical 
z curve which best fit the experimental data was considered as the 
of the matrix and from this value, h, the convection conductance was 
determined. The biggest disadvantage of this method is in the calcula- 
tion of an infinite series for each point and the number of terms neces- 
sary for a good approximation increases as increases. 

In 1950, G. L. Locke (1) developed another method for the evaluation 
of experimental from the theoretical curves. In this method, only 

the maximum slopes of the theoretical and experimental curves need to be 
compared. He derived an expression for the slope of the generalized heat- 
ing curve, (t^ - “ t^) vs * ^/ z > as a function of N^. Locke's 
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expression shows that the maximum slope of the generalized heating curve 
is a unique function of N^: 



d (f2 L i ) 

ln_i< 

dTr7iO 



N 



2 

"tu 



- J x (2i jTN tu ‘ ) e" (N tu + 



T ) 



The above implies that the downstream fluid temperature is measured at 

x = L. Therefore, z = N and t, = t co . This method has become known 
* tu f f2 

as "The Maximum Slope Method". 

One of the idealizations made by Schumann and followed by Locke was 
that the thermal conductivity of the solid parallel to the flow is zero. 
This is very difficult to achieve in an experiment with a continuous flow 
channel in the matrix. Although the effect of finite thermal conductiv- 
ity parallel to the flow is negligible at high Reynold's numbers, it can- 
not be ignored at low Reynold's numbers. Recent works by Crestwick (3), 
Mondt (4), and Howard (5) consider the effects of longitudinal thermal 
conduction in the solution of Schumann's transient heating problem. To 
correlate results, they utilize the conduction parameter, A . This 
parameter is the ratio of conduction heat transfer to the change in in- 
ternal thermal energy of the flowing fluid for the same temperature differ- 
ence. When A = 0, the maximum slope is a unique function of N fcu> as put 
forth by Locke, however, when X f 0, the maximum slope is a function of 

N and A . 
tu 

Howard, in his paper (5), included graphs and a table of maximum 
slopes as a function of N fcu and A . He derived these values by finite 
difference calculations on a digital computer and proposed that the values 
of maximum slopes, as presented, are within + 27 a for a given N and A . 
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Ideally then, if the idealizations, initial and boundary conditions 
can be adhered to, all that is required for the experimental evaluation 



of the heat transfer characteristics of a matrix is to: 

1) Determine the conduction parameter A . 

2) Experimentally determine the maximum slope of the heating or 
cooling curve. 
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3. 



OBJECTIVES. 



In the past, there has been a variety of experiments conceived to 
utilize the principles previously discussed, by approaching the idealiza- 
tions, initial and boundary conditions as closely as possible. Perhaps 
the most difficult condition to achieve is; "at time t = 0 the tempera- 
ture of the incoming fluid changes instantaneously". 

Locke (1) and Coppage (9) both employed an air switching device with 
which air could be quickly directed through the matrix or by-passed to 
the atmosphere. By use of this switch they imposed a "step change" in 
temperature on their test matrices. To approach the "step change" in 
temperature, "requires a fast operating fluid switch for introducing 
the hotter (or colder) fluid in order to assure an instantaneous change 
of fluid temperature. However, even with rapid switching, one of the 
major difficulties here is the inevitable thermal capacity of ducting, 
switch device and matrix holding device causes unavoidable changes in the 
entering temperature". (9) 

Solar Aircraft Company recently did some work testing core modules 
for their T-600 engine (10). In their experimental rig, they utilized 
two steady streams of air and a carriage tray that houses the core module. 
The tray was inserted through the side of the test cabinet and was moved 
from a warm air stream to a ambient or cold air stream to simulate a step 
flunction change of temperature. This type of switching was done manu- 
ally. 

Howard, in his work with new compact surfaces (6), employed a "slid- 
ing drawer". This drawer holds the matrix and the downstream thermo- 
couples. The drawer can be inserted or withdrawn from the test section 



through which ambient air flows continuously. The matrix is heated 
when the drawer is in the withdrawn position by an independent tempera- 
ture controlled air stream. 

One of the objectives of this thesis was to construct and test an 
experimental rig that would improve upon those previously used. This 
test apparatus described herein differs from other test facilities 
chiefly in the method used to simulate the initial condition of a step 
change in temperature of the fluid (air). 

The test apparatus constructed has but one air stream, the air is 
heated by a high resistance nichrome wire heater complex, an integral 
part of the apparatus, located just upstream of the test matrix. The 
time constant or response of the heaters is at least an order of magni- 
tude smaller than the response of the thermocouples and recording devices 
used. By turning the heaters on or off, the step change in temperature 
is simulated. The other systems required some manual dexterity on the 
part of the experimenter to begin a data run, whereas this system in- 
volves only turning off or on the heaters. The temperature recording 
device is started prior to the start of a run. The results achieved by 
this system were an improvement over other systems. With the "sliding 
drawer" and "quick changing valve" there is a change in the flow and there 
fore the hydrodynamic as well as the thermal boundary layer must develope. 
In this system, the hydrodynamic boundary layer is already formed, only 
the thermal boundary layer is affected. A complete description of equip- 
ment is presented in APPENDIX B. 

In addition to the construction of this test apparatus and the verifi 
cation of its reliability by comparing with published data, extension of 
the curves published by Howard (6) into the high Reynold's number region 
was attempted. 
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Finally, heat transfer and flow friction characteristics of several 



new surfaces were determined. 
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4. EXPERIMENTAL TECHNIQUES 

In reference (5) Howard presents results from a numerical finite- 
difference method of calculating the transient behavior of a porous 
media when subjected to a step change in fluid temperature. These cal- 
culations consider the case where longitudinal thermal heat conduction 
cannot be neglected. These data are intimately involved in the determina- 
tion of N fcu and subsequently the Colburn j factor as presented in this 
thesis. For this reason the restrictions, idealizations and conditions 
Howard placed on his solution will be restated here. It is these re- 
strictions, idealizations and conditions that the experimental set-up 
must approach as closely as possible to achieve reliable results. 

1) The flow entering the matrix is steady and uniform in velocity 
and temperature, and remains steady and uniform at any cross 
section as it flows through the matrix. 

2) The thermal conductivity of the matrix is infinite in the 
direction normal to the flow and finite in the direction para- 
llel to the flow. 

3) The thermal capacity of the fluid contained at ary time within 
the matrix is small compared with the thermal capacity of the 
matrix. 

4) The thermal properties of the fluid and matrix are constant 
and uniform. 

5) The convective heat transfer coefficient is some suitable 
average and remains uniform and constant. 

6) At time zero the change in fluid temperature will be a step 
function, with the matrix and its entrained fluid initially at 
some uniform and equal temperature. 

The first idealization was attained by careful experimentation. 
Uniform velocity profiles were obtained by wire screen flow straighteners . 
Uniform temperature profiles were obtained by a wire screen type heater 
explained in detail in APPENDIX B. The velocity profiles were taken at 
the inlet of the matrix over entire range of flow rates to ensure that 



11 



the profiles were uniform. The temperature profiles were taken over the 
entire range of flow rates and found to be uniform within + .5°F. 

The second idealization cannot be strictly realized in practice. The 
effect of the deviations from this idealization can be made negligible in 
the solid if the resistance to heat transfer within the solid is small 
compared to the heat transfer resistance between the solid and fluid. 

The third idealization restricts the fluid to a gas rather than a 
liquid. 

The fourth idealization can be realized if the temperature variation 
is limited to 20 - 25°F. rise above ambient. This would place the 
operating temperature range at 60 to 80°F. In this range the thermal 
properties of the fluid, particularly a gas, vary insignificantly and can 
be considered, as a good approximation, constant and uniform. 

The fifth idealization is imposed upon the computer solution and, of 
course, the determination of the convective heat transfer coefficient is 
the ultimate result of the experimental process. 

Finally, two sets of thermocouples were placed one upstream and one 
downstream of the matrix. Prior to the beginning of each run, these 
thermocouple complexes were read. When a steady temperature differential 
was indicated, the run was begun. 

Herein lies the improvement over the "sliding drawer" technique. A 
series of nichrome wire heaters were installed in the air stream to heat 
the air and subsequently the matrix. To begin a run, the heaters were 
turned on or off, and the temperature vs. time recorder started. The time 
constant (TC) of the nichrome wire heater is very small (TC = .0425 sec 
for m = 1000 lb/hr, TC = .259 sec for m = 10 lb/hr) and hence the approach 



12 



to the final idealization: at time t = 0 the change in fluid temperature 

will be a step function. 

A schematic diagram of the experimental system in which air is used 
as the working fluid is shown in Figure 1. Figure 2 is a photograph of 
the entire test apparatus. Essentially, the system consists of a source 
of air, an air metering device, flow straightener screen pack, an air 
heater section, a test matrix, a pressure measuring system and a tempera- 
ture measuring system. 

Air is drawn through the apparatus by a six-stage centrifugal air com- 
pressor. Air flow rates are measured by an ASME standard D and h/2 ori- 
fice meter with changeable orifice plates. Mass flow through the system 
is controlled by two sliding valces, one on the inlet to the compressor 
and one on the compressor outlet. In addition, there is a fine control 
by-pass valve located just ahead of the compressor. As the air passes 
through the orifice plate it comes to the screen flow straightener at 
which point the flow cross-section changes from round to square. The air 
then flows through the air heater section which serves a dual purpose of 
straightening the air velocity profile and heating the air. In this 
section the air is heated uniformly to a predetermined temperature, usually 
20°F above ambient. The heated air passes into the matrix and heats the 
matrix to the same uniform temperature. At this point, a run is ready to 
start. By referencing the thermocouples downstream of the matrix to the 
thermocouples upstream of the heaters, the initial temperature differen- 
tial can be very closely controlled. This temperature differential is 
continuously recorded by a Leeds and Northrup "Speedomax" strip chart 
recorder when a run is started. The recorded trace has a distinct ad- 
vantage in that there is no transposition of data necessary, hence, a 
decrease in the uncertainty of the results. The resultant trace or curve 



13 



is the primary data for determining j. The N^ u of the matrix can be evalu- 
ated by determining the maximum slope of this fluid temperature differ- 
ence curve during the cooling transient. This is the basis of the maxi- 
mum slope technique as developed by Locke (1). 

The flow friction performances of the matrices tested were obtained 
by determining the isothermal pressure drop across the matrix and the 
mass flow rate of the fluid through the matrix. The pressure drop data 
were obtained from static pressure taps located in the test section im- , 
mediately upstream and downstream of the test section. The entering 
fluid static pressure was also recorded. The mass flow rate data were 
obtained in the air metering section as previously described. In summary, 
the following data were recorded for each run: 

a) Orifice diameter (D ) 

b) Pressure drop across the orifice ( AP ) 

c) Pressure drop across the matrix ( AP ? 

d) Entering static pressure (Ps) 

e) Orifice temperature (T ) 

f) Atmosphere pressure (P ) 

g) Temperature-Time cooling or heating curve 

h) Temperature recorder chart speed 

A complete description of the experimental apparatus and equipment is 
in APPENDIX B. APPENDIX A discusses the relation of the data to the 
equations governing j & f. 
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5. DESCRIPTION OF SURFACES. 

There are data presented in this thesis for seven surfaces. Of 
the seven, three are new. Three of the remaining four surfaces were 
tested by Howard (6). The remaining surface was used primarily as a 
test surface to check the reliability of the experimental apparatus. 

The test surface was a matrix consisting of 20 randomly stacked 
stainless steel, 60 x 6C mesh screens. A photograph of the matrix 
along with a tabulation of its physical and geometrical properties may 
be found in Figure 3. The heat transfer and flow friction characteris- 
tics for such a matrix were published by Coppage in 1952 9), and was 

considered a good check on the performance of the experimental set-up. 

The three matrices examined by Howard and re-examined in this thesis 

were : 

1) Plate-fin reference matrix 

2) 10° skew passage matrix 

3) 20° skew passage matrix 

These particular surfaces were chosen to provide an additional check 
on the reliability of the experimental setup and to provide a compari- 
son between the "integral heater" and "sliding drawer" experimental techni- 
ques. 

The plate-fin surface is formed by corrugated sheets separated by a 
flat spacing sheet to prevent meshing of the corrugations. This matrix 
was to serve as a reference, since the theoretical laminar flow solutions 
for both heat transfer and frictional behavior of triangular cross-sec- 
tions are available. However, difficulty was encountered in forming the 
material into triangular cross-sections, and therefore it fails as a re- 
ference, but the results were interesting and will be discussed later. 
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Figure 4 presents photographs of the plate-fin matrix and a tabulation 
of its' physical and geometrical properties. 

The skewed passage also employs corrugated sheets, but unlike the 
plate-fin, the spacer plate is not included. To prevent meshing of the 
corrugations, the sheets are skewed relative to one another forming what 
is defined as the "skew angle". A photograph of the matrix and a tabula- 
tion of its properties may be found in Figure 5. 

Of the three new surfaces tested, one is a version of the previous- 
ly tested 20° skew matrix using perforated or slotted nickel plates for 
a material vice a continuous material. The plates were corrugated on 
the same device used for the original 20° skew to maintain the same corru- 
gation geometry. The perforations should reduce stagnant areas in the 
matrix, thereby increasing the heat transfer characteristics. 

Finally, two matrices of a parallel plate design were constructed 
and tested. The "parallel plate" matrix is similar to the plate-fin 
in that a spacer plate is placed between the formed sheets to prevent mesh- 
ing. The sheets are formed by rolling. When assembled, the matrix is a 
maze of rectangular passages with an aspect ratio of 77. Aspect ratio is 
defined in reference (8) as the ratio of the two sides of a rectnagular 
flow passage, b/a. The sketch below is, ideally, what was intended in 
designing the "parallel plate" matrix. 




One matrix was formed from a brass material and is used as a refer- 
ence to theoretical laminar flow data for parallel configurations (8). 

The second matrix was formed from the perforated nickel plate to examine 
the effects of the perforations on this type of geometry. Pictures and a 
tabulation of properties for the "parallel plate" matrices are in Figure 6. 
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6 . 



PRESENTATION OF RESULTS. 



The heat transfer and flow friction data for each matrix are pre- 
sented in Tables I to VII. These data are also presented in graphic 

2/3 

form, the Colburn heat transfer factor, j ( = )> anc * t ^ ie friction 

factor, f, are plotted as functions of Reynold's Number, N rS in Figures 7 
through 16. Plotted also in Figures 13 through 15, for comparison pur- 
poses, are the theoretical laminar flow solutions for parallel plates. 

For the heat transfer data, Prandtl number (N pR ) was not really a 

test variable, since the working fluid was air at moderate temperatures. 
2/3 

N pR is included as an approximation to the effect of Prandtl number 

over a moderate range of Prandtl numbers. Basic laminar boundary layer 
solutions indicate that in the Prandtl number range of 0.5 to 1.0 (gases), 
the Prandtl number enters the solution as approximately the 2/3 power. 

Since a moderate temperature range is used in the experiment, air 
properties throughout the matrix are practically constant. Therefore, 
all the properties of air used in the calculations of the heat transfer 
and friction characteristics were those of air at the bulk temperature of 
the air stream as measured at the orifice inlet. 

The Reynold's number is based on the hydraulic diameter and mass flow 
velocity. 



M- 



The hydraulic diameter is defined as: 



D 



H 



= 4 
= 4 



Flow Passage Cross-sectional area 
Flow Passage Wetted Perimeter 

Void Volume x L 

L Heat Transfer Area 



L 



17 



The mass flow velocity is based on the matrix free flow cross-sectional 
area : 



G = 



pA 



The friction factor includes the entrance, exit and flow accelera- 
tion effects. The entrance and exit effects normally provide only a 
small contribution to the overall pressure drop in a matrix. Core fric- 
tion is the term that controls the magnitude of the pressure drop. As a 
result, high accuracy in the evaluation of the entrance and exit loss co- 
efficients is not required. These coefficients k c and are functions 
of flow geometry and Reynold's number and were obtained from Chapter IV, 
Figures 18 to 21, as appropriate, reference (8). 
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7. UNCERTAINTIES AND LIMITATIONS. 

As put forth in the theory, the experimental method utilized, re- 
quires great care to approach the idealizations and initial and boundary 
conditions of the theoretical analysis. The magnitude of errors intro- 
duced if the experimental system does not duplicate the ideal system is 
difficult to evaluate. Each of the conditions and idealizations were 
carefully considered and the experimental system was designed to approach 
as closely as possible the ideal system. These items are discussed in 
the "Experimental Technique" section. No attempt will be made here to 
evaluate or estimate this type of error. 

Other experimental errors are caused by: 1) uncertainty in the 
value of physical constants, 2) inaccuracy in geometric measurements 
and 3) instrument inaccuracies. The method employed herein to describe 
the uncertainty in the values reported is explained fully in Reference (14). 
However, before stating these uncertainties, an example of the process 
follows: 

A theorem of statistics as stated in Reference (14) is: 

"If R is a linear function of n dependent variables, each of which 
is normally distributed, then the relation between the interval for the 
variables, w^, and the interval for the result, w R , which gives the same 
odds for each of the variables and for the result is: 




where w^ = uncertainty interval based on certain odds. 
Consider the Colburn j factor 



j = N 



A£ n 2/3 

tu A PR 



by the theorem and letting P = N_ 



2/3 
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(Ap %, a ) 2 + w v f + (^tj,4 P w A )*+ (|i^4 p' 1 ^ Wp) 1 



if both sides are divided by j then: 



= 

J 



+ (W 



( 1 ) 



Before an equation such as this can be evaluated the uncertainty 
terms must be known. 

The physical constants necessary are listed below. References (12) 
and (13) were sources for the values of these constants. The uncertainties 
in these values, as closely as can be determined, are also listed below. 

c + 0.5% 

s — 

c + 0.5% 






+ 1 . 0 % 
+ 2 . 0 % 



k + 0.5% 

s - 



The uncertainty in physical linear measurements could be consider- 
ed negligible if core construction did not have minor inconsistencies. 

The uncertainty in the weight of a core is considered negligible. Conser- 
vative estimates of the uncertainty in linear measurements are: 

A, A. , A , A , + 1.0% 
fr* c* s’ - 

L + 0.5% 

W < 0.1% 
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Instrumentation errors were such that only the manometers need be 
considered. With the exception of the orifice temperature, all tempera- 
tures measured were temperature differences and these were measured in 
inches as indicated in Appendix A. Since pressures varied and manometers 
to measure them changed, the maximum uncertainty encountered was used in 
the error analysis and is listed below. 

P + 1.25% 

o — 

P + 1.70% 

m — 

P at r. negligible 

Another uncertainty which requires special mention is that in N^. 

Since j does not vary linearly with N R , the uncertainty at high and low 

N R was considered. This uncertainty was determined by utilizing Figure 19 
and the uncertainty in Maximum Slope as 2.0%. 

IT = 3.5, X= 0 dW tu = + 6. OX 

tu M— " 

tu 

N = 30, A= .04 dN tu = + 9.6% 

tu * — — — 



Returning to Equation (1) and substituting as appropriate: 



For high N : - j 

R J 1 


+ C-3 4 ) 2 


]\ ..13 


tA/j' r 6.3% 
J 


For low N_: = 1 

R J 1 


It + 0 °')' +(• O')" 1 + 


Vz 

^ .098 


wj = %8 % 
J 



Following the above type of analysis, the following additional 
uncertainties were calculated. 

m +1.0% 

N r + 2 . 5% 

f + 4.3% 
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8. DISCUSSION OF RESULTS. 



The 60 x 60 mesh screen matrix was used solely as a test matrix, a 
matrix whose heat transfer and friction characteristics were well docu- 
mented. Figure 7 and Table I present the results. There is quite good 
agreement in the friction factor but a slight variation in j factor. This 
difference is attributed to the difference in the geometries of the test 
matrices. Exact reproduction of a screen matrix is difficult to achieve. 
The correlation was such that the experimental test setup was considered 
capable of delivering reliable data. 

For further verification on the reliability of the test setup. Figures 
8 to 10 and Tables II to IV present the test results for the Plate-Fin 
Reference Matrix, the 10° Skew Stainless Steel Matrix and the 20° Skew 
Stainless Steel Matrix used by Howard (6). Although this was not the 
original purpose for utilizing these matrices, the results do indicate 
conclusively that the "integral heater" is as good or better than the 
"sliding drawer" technique. 

The original purpose for utilizing these matrices was to try and ex- 
tend Howard's curves into the higher Reynold's number region (lower N fcu 
region). The points obtained in this region are plotted on the curves and 
the perinent data are recorded in the Tables. None of them follow the 
straight line of their predecessors. There are at least two reasons for 
this obvious departure from the norm. First, it is a well known fact 
that the "Maximum Slope" Method has an inherently large error in the low 
Nj_ u region of 1. 5<N tu <3. 5. This fact was originally published by Locke 
(1) for A= 0. Howard (5) took into account longitudinal conduction and 
published an error curve that includes the effect of longitudinal conduc- 
tion on an error in N^, see Figure 19. In the regions below N tu <1.5 
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the error is acceptable, yet, a discrepancy still exists. This points 
out the second reason for the discrepancy. The maximum slope at high 
flow rates (low N^) occurs at time t = 0+. Instrumentation prevents 
accurate measurement of the true maximum slope at this time due to in- 
herent delay or lag time. 

Extrapolation of the time- temperature curves back to time t = 0 was 
attempted in an effort to remove the effects of instrumentation. The 
method used was either incorrect or unreliable. Satisfactory results 
were unobtainable. 

The merits of the plate-fin, 10° skew and 20° skew matrices as heat 
transfer surfaces was discussed thoroughly by Howard in reference (6). 

Figures 11, 13 and 14 and Tables V, VI and VII present the results 
for a new geometry tested and a new idea in material. Figure 11 and 
Table V present the results for a 20° skew nickel matrix constructed from 
slotted plates. Appendix C explains fully the geometry of the slots and 
the calculations that determine the basic geometric parameters p, j3 , 
and r^. The effects of longitudinal conduction was considered negligible 
as the matrix was constructed so that the air flow was normal to the long 
dimension of the slots. Plotting all the points uncorrected for conduc- 
tion showed an erratic departure from the expected results for < 200. 
It was then decided that perhaps the basic assumption was not warranted 
completely. Corrections for longitudinal conduction were made where 
significant and these points were plotted. These corrected values of j 
are overestimations, for, although A g is corrected for the effects of the 
slots, the matrix length L is not the true conduction path length but at 
best a lesser approximation. Since X is obviously not equal to zero, 
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the true value of j is probably between the uncorrected j ( A = 0) and 
the j corrected for longitudinal conduction using L as an approximation 
for conduction path length. This was considered in the interpretation of 
results . 

Figure 12 compares the 20° skew stainless steel matrix with the per- 
forated nickel matrix. There is a 70 to 100% increase in j, while f in- 
creases 10 to 25%. 

It must be understood that these two matrices are geometrically similar 
only in overall dimensions. Actually, the heat transfer areas were cal- 
culated on different bases. The 20° skew stainless steel heat transfer 
area was based on the total surface area available for heat transfer while 
the heat transfer area of the 20° skew perforated was based on this same 
area less the area of the slots. Also, not included in the heat trans- 
fer area is the peripheral area of the slots normal to the flow. If 
the 20° skew perforated j factor were based on the total available area 
for heat transfer, the values plotted would be decreased approximately 
40%. The perforated material still represents a substantial improvement 
in the performance. This improvement is attributed to the reduction of 
inefficient stagnant areas within the matrix and the alteration of the 
thermal and hydrodynamic boundary layers caused by the slots. 

Figure 13 and Table VI present the results for a "parallel plate" mat- 
rix constructed from brass shim stock. As can be seen in Figure 13, the 
heat transfer performance falls far below that of theoretical laminar 
flow performance for parallel plates. Several things may account for this: 
the aspect ratio of the flow channels in this matrix was 77 vice co for 
truly parallel plates, and also, the flow channels were not uniform in 
cross-section from one to another. The non-uniformity in cross-section 
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results in uneven mass flow through the matrix, hence inefficient use 
of the available heat transfer surfaces, and therefore a decrease in j. 

The decrease in friction factor is possibly due to the sharp angle 
regions where the formed sections contact the flat plate. The wall shear 
stress in this type of region is small which would tend to reduce the 
value of f. These regions would also tend to produce stagnant regions 
previously mentioned, having a diminishing effect also on j. 

Figure 14 and Table VII present the results for a "parallel plate" 
matrix constructed from perforated nickel plate. This matrix was construct- 
ed in the identical fashion as the brass matrix. Although it is not read- 
ily apparent in Figure 14, the heat transfer characteristics have been 
improved nearly 80%, while the friction factor has remained the same. 

Again it is necessary to make note of the fact that the heat transfer area 
of the two matrices were calculated on different bases. These bases are 
the same as discussed for the 20° skew matrices. Applicable also is the 
approximate reduction by 40% of the plotted values of the j factor if 
all of the area, in the perforated material, available for heat transfer, 
were considered in the calculations. 

Figure 15 is a comparison of the two "parallel plate" matrices. The 
improvement realized by using perforated material is readily apparent, 
even considering the 40% reduction mentioned above. The additional per- 
formance gain in the transferral of heat is again attributed to the re- 
lieving of the stagnant areas and the alteration of the boundary layers 
caused by the slots. 

The fact that the friction factor did not increase does not lend 
itself to previous conclusions (20° skew) ; that the presence of slots 
increases the friction factoi . Some additional factor must be necessary. 
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Perhaps, the .skewing coupled with the slots contribute more to friction 
then straight, slotted channels. 

In the determination of the j factor for the perforated "parallel 
plate" matrix the conduction parameter, A , was considered zero. Here 
again the plotted points departed from the expected straight line indicat- 
ing that X , in fact, was not zero. Corrections for longitudinal conduc- 
tion were made using the matrix length, L, as an approximation for the 
conduction path length. These points were plotted also. Since the cor- 
rected points are overestimations and A / 0, the true value of j lies 
between these limits. 

The designer of a compact heat exchanger is generally restricted as to 
frontal area, volume, weight and pressure drop. There is presently coming 
into use, several means of comparing heat transfer surfaces with these 
restrictions in mind. Essentially, various parameters are compared 
graphically so that the relative merit of the various surfaces may be es- 
tablished. 

The ratio j/f is plotted vs. Reynold's Number as a means of relating 
frontal flow area with pressure drop. For a given pressure drop a high 
j/f characteristic implies a small flow frontal area, j/f when plotted 
vs. N r is known as the Flow Frontal Area "Goodness" factor. 

The heat transfer power, h gTD> is plotted vs. the flow friction power, 
E STD’ 3n< * t * ie resu l tant characteristic curve is a measure of relative core 
volumes. Generally, the surface with a high h vs. E characteristic 

□iD ui D 

will tend to have a small core volume, h ^ plotted in this fashion is 
known as the Volume "Goodness" factor. 

Figure 16 compares the 20° skew matrices and the "parallel plate" 
matrices. Although the "Goodness" factors are not infallible, they are 
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indicative of relative performance. Tn this light, the perforated 
material should be excellent matrix material. 
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9. 



CONCLUSIONS. 



Based upon the correlation of present data with that obtained by 
previous experimenters, the experimental test rig design presented here- 
in is considered capable of producing reliable test data on the heat trans- 
fer and friction characteristics of a core module. 

Basic heat transfer and flow friction characteristics for a new 
geometry and a new idea in material are presented in Figures 11 to 15. 

The geometric and physical properties are summarized in Figures 5 and 6. 
There appears to be a definite advantage to constructing matrices from 
perforated material. Performance is improved considerably and based on 
the "Goodness" factors (Figure 16) size and weight advantages are realiz- 
ed. 

The 20° skew perforated nickel matrix was the top surface in perfor- 
mance. This matrix had the highest Flow Frontal Area "Goodness" factor, 
as well as the highest Volume to surface "Gcodness" factor. This does not 
mean that this is the "optimum" surface for any given application but it 
does merit consideration. 

In both geometries tested (20° skew and "parallel plate") the per- 
forated material out-performed its counterpart. This would also indicate 
that further testing of this type of material (slotted) is warranted. 

It would seem that this material and some geometry should result in an out- 
standing heat exchanger surface. 
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10. RECOMMENDATIONS FOR FUTURE WORK. 



1) Some means should be found to reliably determine N tu in the 

high N r range (1.0 3.5). Most promising is the use 

of a periodic or cycling technique in which effectiveness is 

measured and used to find N rather than the maximum slope. 

tu r 

This ability is necessary if performance curves are to be extend- 
ed to N r “>■ 500. 

2) More work can be done with perforated material. Different 
geometries should be tested to see if consistantly high per- 
formances are realized as they have been in the past. Addition- 
ally s this type of surface is suspect of fouling. If these 
materials are as good as they seem to be, then examination of 
the effects of fouling on their performance is in order. 

3) As more surfaces become available, they should be tested 

and their heat transfer and friction characteristics determined. 
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FIGURE 1. SCHEMATIC DIAGRAM OF TEST APPARATUS 




18-8 St. 

Matrix Material Steel 

specific Heat (o ) BTU/lh °F °- 12 

Thermal Conductivity (k s > BTU/hr °F ft 10.0 



Material Thickness, inches ^ 

Total Heat Transfer Area (A) ft 
Frontal Area (A fr ) ft 2 

Total Conduction Area (A g ) ft 
Flow Cross Section Area (A c ) ft 
Matrix Volume (V m ) ft ^ 
Matrix Density lh/ft 



0.017 

4.0159 

0.07055 

0.02215 

0.0484 

0.00199 

155. 



Hydraulic Diameter (D R ) ft 
Compactness (y0) ft^/ft^ 
Porosity (p) 



0.00156 

2009. 

0.686 




FIGURE 3 



Physical and Geometric Properties of 
60x60 Mesh Screen Matrix 
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Type 302 



Matrix Material 


S. Steel 


Specific Heat (c s ) BTU/lb °F 


0.11 


Thermal Conductivity (k g ) BTU/hr °P ft 


7.0 


Material Thickness, inches 


0.001 


Total Heat Transfer Area (A) ft 2 


30.12 


Frontal Area (A fr ) ft 2 


0.07023 


Total Conduction Area (A Q ) ft 2 


0.00752 


Flow Cross Section Area (A c ) ft 2 


0.06271 


Matrix Volume (V m ) ft^ 


0.01170 


Matrix Density (p m ) lb/ft^ 


52.9 


Hydraulic Diameter (D^) ft 


.001388 


Compactness (j3 ) ft 2 /ft^ 


2573* 


Porosity (p) 


0.893 



FIGURE 4 Physical and Geometric Properties of 
Plate-Pin Reference Matrix 
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Matrix Material 

Specific Heat (d g ) BTU/lb °F 

Thermal Cond. (kg) BTU/hr °F ft 

Material Thickness, inches 

2 

Total Heat Transfer Area (A) ft 

Frontal Area (A-,) ft 2 

fr o 

Total Conduction Area (A e ) ft 

s o 

Flow Cross Section Area (A c > ft 
Matrix Volume (V m ) ft 5 
Matrix Density (|0 m ) lh/ft^ 

Hydraulic Diameter (Dg) ft 
Compactness (/3) ft 2 /ft^ 

Porosity (p) 



10° SKEW 20° SKEW 20° SKEW 



Type 302 
S. Steel 


Type 302 
S. Steel 


perf . 
Nickel 


0.11- 


0.11 


0.1065 


7.0 


7.0 


38.7 


0.001 


0.001 


0.0022 


18.38 


18.04 


15.41 


0.07024 


0.0702 


0.0734 


0.00464 


0.00452 


0.00349 


0.0656 


0.0657 


0.06412 


0.0117 


0.01157 


0.01211 


32.4 


32.22 


50.1 


.002381 


.002416 


.002528 


1570. 


1549. 


1272. 


0.934 


0.936 


0.873 



FIGURE 5 Physical and Geometric Properties of 
Skew Matrices 
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Matrix Material 


70-30 

Brass 


Perf. 

Nickel 


Specific Heat (c s ) BTU/lb °F 


0.092 


0.1065 


Thermal Conductivity Ck g ) BTU/hr °F ft 


66. 


38.7 


Material Thickness, inches 


0.001 


0.0022 


Total Heat Transfer Area (A) ft 2 


22.057 


18.586 


Frontal Area (A fl> ) ft 2 


0.06986 


0.0695 


Total Conduction Area (A g ) ft 2 


0.00551 


0.00416 


Flow Cross Section Area (A Q ) ft 2 


0.06434 


0.0584 


Matrix Volume (V m ) ft 5 


0.01164 


0.01158 


Matrix Density (^ m ) lb/ft^ 


41.5 


63.6 


Hydraulic Diameter (D H ) ft 


.001944 


.001928 


Compactness (0) ft 2 /ft^ 


1895. 


1605. 


Porosity (p) 


.921 


.840 



FIGURE 6 Physical and Geometric Properties of 
"Parallel Plate" Matrices 
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Figure 9 10° Skew Stainless Steel Matrix 

I 
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Figure 10 20° Skew Stainless Steel Matrix 

) 
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Figure 11 20° Skew Nickel (Perforated) Matrix 

) . ' : 
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Figure 12 Comparison of perforated 

and non-perforated 20° Skew 
. Matrices 



4-3 




Figure 13 ."Parallel Plate" Brass Llatrix 
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Figure 14 "Parallel Plate" Nickel (perforated) Matrix 

) . 
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Figure 15 Comparison of Perforated and 

non-perf orated "Parallel Plate" 
. . Matrices 
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200 



0,5 



0.3 



o.i LIJ 2‘<y ^“^o- lim 56 rob' ‘ ,J 200 



(7)20° Skew Perforated Nickel 
(z) "Parallel Plate" Perforated Nickel 
(5) 20° Skew Stainless Steel 
( 4 ) "Parallel Plate" Brass 



Figure 16 Design "Goodness" Factors 
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FIGURE 17 GENERALIZED COOLING CURVE 

NOTE: All measurements are made In inches 
to facilitate computations. 
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NUMBER OF TRANSFER UNITS VERSUS MAXIMUM SLOPE 
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Figure 18 Number of Transfer Units (N. 
Maximum Slope 




S3 





CT> 

H 



50 



Maximum Slope vs 




FI3UEZ 20. kstrix EclEer "psirssn tie* 




-**-rs-rZ 22. lLc~rix nclaer i.r<rrr,5tre£r ;rr 





FIGURE 22. Matrix Holder modified with probe 
for velocity profile measurement 
and temperature distribution 
measurement (Downstream View) 





FIGURE 23. Nichrome Wire Heater (Upstream View) 




FIGURE 2lu Nichrome Wire Heater (Downstream View) 




PROGRAM MAXSLO 
L =2 
K= 1 

30 READ 300»XML»ASfRHfP0RfAFLfAHT f WM 
300 FORMAT ( 7F 1 0 • 5 ) 

40 READ 400,EK,CK,SK,CM, N 
400 FORMAK4F10.5, 12) 

M=1 

10 READ 100,DO,DELPO,GAM,ATMP,COR, UH.UFR 

100 FORMAK7F10.5) 

READ 101,CS,DELH , HS , SLO , TEMPO 

101 FORMAT ( 5F 1 0 . 5 ) 

EMDOT=359.0*DO*DO*SQRTF(DELPO*GAM*ATMP/29.92)*COR 
C P = 0 • 24 

RNUMP=S.09*EMDOT/UFR 
RNUMH=J4.0*XML*EMD0T)/( AHT*UH) 

RNUMFR=RNUMH*UH/UFR 

CAPFL=EMDOT*CP/3600.0 

CAPS=WM«CM 

SL0MAX=(CAPS*SL01/ (CAPFL*CS1 
CONDPAR=(SK*AS)/( XML*EMC0T*CP ) 

B=CK+EK 

C= 1 . 0 + P0R» «2 

G=EMD0T/(AFL*3600.0> 

DELP=DELH*. 03613 
PS=HS*. 03613 
PA=ATMP*.4912 
PI =PA-PS 
P2=P 1-DELP 
PM=( P1+P2J/2. 

RHOM= { PM*144.0) / ( 53 . 3 * ( TEMPO+460.0) ) 
FF1=((64.4*RHQM*DELP*144.0)/(G*C) ) 

FF2=C*DELP/PM 



FF3=FF1 -FF2-B 
FFR=FF3*RH/XML 
WRITE OUTPUT TAPE 
FORMAT ( 9H RUN NO. 
WRITE OUTPUT TAPE 



4,4 ,M 
12 /// ) 

4,1»EMD0T, RNU MP » R NUMH * CAP FL» CAPS 

1 FORMAT ( 7H EMDOT= F 1 0. 5 , 3X , 7H RNUMP=E 20. 8 , 3X , 7H RNUMH=F10. 
17H CAPFL=F10.5,3X,6H CAPS=F1 0.5///) 

WRITE OUTPUT TAPE 4, 2 , CONDPAR , S LOMA X 

2 FORMAT ( 9H CONDP AR=F 1 0 . 5 , 3X , 8H S LOMA X=F 1 0. 5/ // ) 

WRITE OUTPUT TAPE 4 , 3 , G , RNUMFR f FFR 

3 FORMAT ( 3H G=F 1 0 .5, 3X , 8H RNUMFR=F10. 5, 3X, 5H FFR=F10.5/1 
M=M + 1 

IF(N-M)50, 10,10 
50 K=K* 1 

I F ( L— K 120,30,30 
20 STOP 
END 

END _ . 



Figure 27 Program MaxSIo Readout 
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TABLE I SUMMARY OP HEAT TRANSFER AND FRICTIONAL RESULTS 
60x60 MESH STAINLESS STEEL SCREEN MATRIX 



HEAT TRANSFER RESTTT.TR . 

RUN REYNOLD’S N Q1TI N TO 
/ NUMBER 

1 431.57 

2 234.81 

3 159.26 

4 129.69 

5 106.78 

6 79.88 

7 50.69 

8 29.60 

9 15.43 

10 8.49 



CONDUCTION 


MAXIMUM 


N , 


PARAMETER 


SLOPE 


tu 


0.0 


.58333 


3.12 


0.0 


.68566 


4.96 


0.0 


.74639 


6.09 


0.0 


.80395 


7.25 


0.0 


.82926 


7.77 


0.0 


.87817 


8.83 


0.0 


1.05579 


11.90 


0.0 


1.22105 


17.94 


0.0 


1.57029 


30.19 


0.0 


2.11897 


55.65 



.0267 

.0425 

.0521 

.0620 

.0659 

.0755 

.1040 

.1535 

.2584 

.4760 



ISOTHERMAL FRICTTON 



RUN REYNOLD’S 
# NUMBER 

1 437.93 

2 238.27 

3 161.60 

4 131.60 

5 108.36 

6 81.07 

7 51.44 

8 50.04 

9 15.66 

10 8.62 



RESULTS * 

FRICTION 

FACTOR 

.68548 

.82196 

.89263 

.97400 

1.06410 

1.23243 

1.57979 

2.31547 

4.47090 

9.50629 
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TABLE II SUMMARY OP HEAT TRANSFER AND FRICTIONAL RESULTS: 
PLATE-FIN REFERENCE MATRIX 



HEAT 


TRANSFER 


RESULTS: 










RUN 

# 


REYNOLD'S 

NUMBER' 

% 


N N 
W ST W PR 

5 


CONDUCTION 

PARAMETER 

X 


MAXIMUM 

SLOPE 


N tu 


• N tu 

U=o) 


1 


505.17 


.00277 


.00128 


.52585 


1.68 


1.68 


2 


433.26 


.00429 


.00150 


.55712 


2.60 


2.60 


3 


387.19 


.00505 


.00167 


.58044 


5.06 


5.06 


4 


522.08 


.00654 


.00201 


.62585 


3.84 


■ 3.84 


5 


248.92 


.00842 


.00260 


.69377 


5.10 


5.10 


6 


187.25 


.01162 


.00546 


.79434 


7.04 


7.04 


7 


125.74 


.01761 


.00515 


.94678 


10.67 


10.44 


8 


75.04 


.02889 


.00862 


1.16341 


17.51 


16.21 


9 


42.68 


.05699 


.01516 


1.43365 


34.54 


25.07 


10 


25.81 


.08405 


.02717 


1.43851 


50.93 


25.22 



ISOTHERMAL FRICTION RESULTS: 



RUN 

# 


REYNOLD'S 

NUMBER 

n r 


FRICTION 

FACTOR 

f 


1 


512.59 


.02197 


2 


439.62 


.02485 


3 


392.86 


.02795 


4 


326.79 


.05304 


5 


252.57 


.04210 


6 


189.99 


.05433 


7 


127.58 


.08205 


8 


76.14 


.12974 


9 


43.30 


.25164 


10 


24.17 


.39657 
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TABLE III SUMMARY OP HEAT TRANSFER AND FRICTIONAL RESULTS: 
10° SKEW STAINLESS STEEL MATRIX 



HEAT 


TRANSFER 


RESULTS: 










RUN 


REYNOLD’S 


N N 
oT rR 


CONDUCTION 


MAXIMUM 


N tu 


N tu 


# 


NUMBER 


PARAMETER 


SLOPE 


% 


3 


A 






(A =0) 


1 


805.27 


.00375 ‘ 


.00081 


.46533 


1.32 


1.32 


2 


643.52 


.00440 


.00101 


.51088 


1.55 


1.55 


5 


517.86 


.00695 


.00126 


^55129 


2.45 


2.45. 


4 


402.58 


.01073 


.00162 


.62114 


3.78 


3.78 


5 


305.00 


.01488 


.00214 


.70167 


5.24 


5.24 


6 


200.72 


.02271 


.00324 


.83894 


8.00 


7.98 


7 


119.58 


.03642 


.00544 


1.02840 


12.83 


12.46 


8 


68.74 


.06376 


.00953 


1.27448 


22.46 


19.62 


9 


37.61 


.10788 


.01728 


1.44693 


38.00 


25.53 



ISOTHERMAL FRICTION FACTOR: 



RUN 

* 


REYNOLD’S 

NUMBER 

n r 


FRICTION 

FACTOR 

f 


1 


817.04 


.02791 


2 


652.71 


.03405 


3 


525.43 


.04092 


4 


408.46 


.04951 


5 


309.46 


.06345 


6 


203.79 


.09412 


7 


121.41 


.15138 


8 


70.00 


.27472 


9 


38.16 


.47276 
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TABLE IV SUMMARY OP HEAT TRANSFER AND FRICTIONAL RESULTS: 
20° SKEW STAINLESS STEEL MATRIX 



HEAT 


TRANSFER 


RESULTS: . 










RUN 

# 


REYNOLD’S 

NUMBER 

n r 


N N 
"srpR 

j 


CONDUCTION 

PARAMETER 

A 


MAXIMUM 

SLOPE 


N tu 


*tu 

(A=0) 


1 


835.38 


.00394 


.00078 


.47359 


1.36 


1.36 


2 


728.65 


.00414 


.00090 


.48768 


1.43 


1.43 


3 


636.57 


.00475 


.00103 


.52394 


1.64 


1.64 


4 


535.10 


.00805 


.00122 


.56474 


2.78 


2.78 


5 


412.85 


.01182 


.00158 


.63657 


4.09 


4.09 


6 


312.06 


.01524 


.00209 


.70256 


5.26 


5.26 


T 


208.65 


.02167 


.00313 


.81581 


7.48 


7.48 


8 


126.47 


.03676 


.00517 


1.02339 


12.69 


12.33 


9 


71.21 


.06593 


.00917 


1.28940 


22.76 


20.09 


10 


39.79 


.09441 


.01646 


1.39790 


32.59 


23.79 



ISOTHERMAL FRICTION FACTOR: 



RUN 

/ 


REYNOLD’S 

NUMBER 

n r 


FRICTION. 

FACTOR 

f 


1 


848.11 


.03197 


2 


739.39 


.03471 


3 


645.95 


.03908 


4 


542.97 


.04465 


5 


418,94 


.05431 


6 


316.65 


.06856 


7 


211.73 


.10255 


8 


128.33 


.16015 


9 


72.26 


.28784 


10 


40.40 


. .48428 



TABLE V SUMMARY OP HEAT TRANSFER AM) FRICTIONAL RESULTS 
20° SKEW NICKEL MATRIX (PERFORATED PLATES) 



HEAT TRANSFER RESTTT/rs , 



RUN REYNOLD'S 

0 NUMBER 

^r 

1 966.07 

2 835.43 

3 738.42 

4 619.79 

5 484.37 

6 362.97 

7 239.79 

8 145.03 

9 82.26 

10 46.22 

11 23.50 



N N ' 

sr'pR 


CONDUCTION MAXIMUM 


N. 


N 


J 


PARAMETER 

X 


SLOPE 


tu 


tu 

(X=0) 


.00900 


.00337 


.56394 


2.72 


2.72 


.01146 


.00390 


.60157 


3.42 


3.46 


.01219 


.00441 


, .61776 • 


3.70 


3.68 


.01418 


.00526 


.64762 


4.24 


4.28 


.01802 


.00674 


.71212 


5.41 


5.44 


.02497 


.00898 


.81813 


7.58 


7.54 


.03517 


.01359 


.95653 


11.29 


10.62 


.05199 


.02243 


1.14669 


19.65 


15.70 


.0*7879 


.03951 


1.39808 


_ 


23.79 


.08147 


.07039 


1.42079 


_ 


24.60 


.09134 


.15843 


1.50225 




27.58 



ISOTHERMAL FRICTION Pamvw . 



RUN 


REYNOLD'S 


PRICTION 


0 


NUMBER 


FACTOR 




Nr 


f 


1 


980.29 


.03316 


2 


847.72 


.03709 


3 


749.28 


.04116 


4 


628.91 


.04686 


5 


491.63 


.05802 


6 


368.31 


.07208 


? 


243.32 


.10258 


8 


147.16 


.16737 


9 


83.46 


.29667 


10 


46.89 


.49889 


11 


23.85 


1.09090 
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TABLE VI SUMMARY OF HEAT TRANSFER AND FRICTIONAL RESULTS’*. 
"PARALLEL PLATE" BRASS MATRIX 



HEAT TRANSFER RESULTS: 



RUN 

i 


REYNOLD’S 

NUMBER 

N R 


N N 
1 ox PR 

J 


CONDUCTION 

PARAMETER 

X 


MAXIMUM 

SLOPE 


N tu 


N tu 

(A=0) 


1 


683.97 


.00376 


.00896 


.49410 


1.44 


1.46 


2 


551.23 


.00527 


.01112 


.54861 


2.27 


2.40 


3 


437.55 


.00748 


.01401 


.59598 


3.22 


3.36 


4 


340.15 


.01038 


.01802 


.66717 


4.47 


4.63 


5 


255.77 


.01335 


.02390 


.73472 


5.75 


5.63 


6 


168.67 


.02148 


.03634 


.87166 


9.25 


8.71 


7 


102.18 


.03483 


.05999 


.97562 


15.00 


11.15 


8 


57.06 


. - 


.10761 


1.07854 


- 


- 


S) 


33.37 


- 


.18369 


.98829 


- 





ISOTHERMAL FRICTION FACTOR : 



RUN 

* 


REYNOLD’S 

NUMBER 

n r 


FRICTION 

FACTOR 

f 


1 


694.19 


.02846 


2 


559.46 


.03454 


3 


444.09 


.04292 


4 


345.24 


.05240 


5 


259.52 


.06907 


6 


171.19 


.10680 


7 


103.71 


.17239 


8 ' 


57.90 


.32370 


9 


33.86 


.53085 
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TABLE VII SUMMARY OP HEAT TRANSFER AND FRICTIONAL RESULTS: 

"PARALLEL PLATE" NICKEL MATRIX (PERFORATED PLATES) 



HEAT 


TRANSFER 


RESULTS: 










RUN 

* 


REYNOLD' 

NUMBER 

*R 


5 N N 

6 « ST « PR 

3 . 


CONDUCTION 

PARAMETER 

A 


MAXIMUM 

SLOPE 


W tu 


N tu 

( A = o ) 


1 


819.49 


.00590 


.00393 


.54781 


2.32 


2.36 


2 


647.21 


.00810 


.00497 


; 58864 


3.18 


3.24 


3 


518.89 


.00988 


.00620 


.62923 


3.89 


3.95 


4 


400.50 


.01458 


.00804 


.73180 


5.77 


5.83 


5 


304.49 


.01891 


.01058 


.81941 


7.61 


7.56 


6 


200.68 


.02659 


.01602 


.95459 


11.31 


10.63 


7 


121.15 


.04039 


.02658 


1.16142 


21.61 


16.15 


8 


67.53 


.05502 


.04768 


1.34786 


- 


22.00 


9 


37.92 


.04945 


.08478 


1.27951 


- 


19.70 



ISOTHERMAL FRICTION FACTOR: 



RUN 

§ 


REYNOLD’S 

NUMBER 

n r 


FRICTION 

FACTOR 

f 


1 


831.56 


.02395 


2 


656.72 


.02955 


3 


526.51 


.03577 


4 


406.49 


.04439 


5 


309.04 


.05682 


6 


203.63 


.08725 


7 


122.96 


.14242 


8 


68.54 


.26776 


9 


38.48 


.4587 f 
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APPENDIX A. 



DATA REDUCTION RELATIONS 

The following summary of some of the data reduction relations is 
given to clarify some of the definitions of the geometrical parameters, 
to show how they influence the deduced results and to facilitate an 
understand as to the exact procedure required to produce from raw data 
a point on the j vs N^ and f vs N^ curves. 

The three most significant geometric parameters that afford a means 
of comparison between matrices are so defined that there exists a unique 
relation between them. Independent measurement of any two establishes 
the third or serves as a check on its independent measurement. The 
three primary surface geometry parameters of a matrix are: 

1) HYDRAULIC DIAMETER: 







% - 4 r L 


(A-l) 


2) 


POROSITY: 


Matrix Flow 

Void Volume A L A 






P 


= Total Volume = = ~ 

rr A fr 


(A-2) 


3) 


AREA COMPACTNESS: 


Total Heat Transfer 
Surface Area 






h 


Total Volume = t — — 

A fr L 


(A-3) 



NOTE: If equation (A-2) is divided by (A-3) the unique relation mentioned 

previously is formed. 

r R = p/0 (A-4) 



MAXIMUM SLOPE : 

As previously discussed, through the work of Locke, Mondt, Crestwick 
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a.d ’-toward, it has been established -:h;- •: the axi i slope of a general- 



ized heating or cooling curve is a unique function of and A 
conduction parameter, 



fc f 2 ” fc i 

d( t : ~r ) 

fl i 



= f(N' * X ) 



T = generalized time variable 



T 



: air temp, downstream 



T/n , 



me /W c © 
p s s 



and d( T /N fcu ) = tnc^/ W g C g (d © ) 



f2 "fl + 1 
C fl ■ l i 



t 7 - t. 

and: d(-^ p— ) 

fl " i 

combining (A-5) and (A-6) 



t„). 



d( ' 



fl i 
d( T/N. „) 



d( ' 



, the 



(A-5) 



(A-6) 



(A- 7) 



66 



The matrix is heated to a uniform temperature, t^. The heat source 
is secured and the difference between downstream and upstream temperatures 
(t^ " tfj) recorded vs. time. The maximum rate of change of this 
difference with time may be determined from the record. 

From Figure 17 : 






1 

x 



X 

Chart speed 



d & sec 



y = d(t f2 - t fl ) 



inches 



M = (t fl - t^ inches 



W c Matrix heat Capacitance 

* S - , S ' = Fluid heat Capacitance 

mc p 



7T BTU/°F 

BTU/°F sec 



Finally, Equation (A-7) becomes: 



d( fc f2 " fc i ) 
fc fl ' t i 
d( T /N tu ) 



max 



M 



Chart speed 



With this value of the maximum slope and a calculated value of A , 
Howard's curves or tables may be entered to obtain a value of N 

tu 



MASS FLOW FATE : 

The mass flow rate is calculated as outlined in reference (7): 
m = 359 K D* FaY | A P Q t 
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where 



K = flow coefficient, velocity of approach factor 



Fa = thermal expansion factor 

Y = expansion factor 

D = orifice diameter, inches 
o ’ 

^ = specific weight of fluid flowing at the orifice 
3 

temperature lb/ ft , corrected to existing 
atmospheric pressure 

A P q = Pressure drop across theorifice , inches 

REYNOLD'S NUMBER : 



By definition, the Reynolds number is 



R JUL 

Where G is the mass flow velocity based on the flow mean cross- 
sectional area, A . 



(A- 8) 




Note, if (A-9) is substituted into (A-8) 



(A- 9) 



N 



R 



4m 

^ A fr 



( 



n L 
P 



) 



(A- 10) 



but from (A-4) r H/p = £ = ) (, 

Therefore, for a given range of mass flows, the Reynolds number range 
for a given matrix is proportional to the inverse of its compactness. 



n r c=C '/($ 

For ease of calculation, recall from (A-3) 



P 



substituting into (A-ll) yields 



4Lm 

A M 



(A-12) 



N 



THE FANNING FRICTION FACTOR : 

The following sketch describes the flow system to be considered, 
i a. t 2 



1 

1 


TZZZZZZZZZZZZZZZZnZZL 




! r 


777777777777777777777, 


\ 






1 

| ^ 


i 


7777777777777777777777, 


1 


i 

i 




i * 



For gas flow heat exchanger application, the pressure changes 
from section 1 to a and from b to 2 are very small, relative to the 
total pressure, and therefore ~ 'Vj and ~ 

The relation for the flow stream pressure drop calculation for 
most heat exchanger cores is (8): 



AP - g* V7 

P, P. 









entrance flow core exit 

effect acceleration friction effect 



(A 



Solving (A-13) for f and recalling from (A-l) A/A c = L/r H also 
v = /O , we find: 




Now the perfect gas law may be applied since testing is done with air 
moderate temperatures and pressures. Equation (A-14) becomes for the 
isothermal case: 



-13) 



-14) 

at 
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f. 






T ( A “ 15 ) 



p = and T, = T„ 

/ m 



, _L RT 

where p = p m ~ ^ olIU A l " "2 

The first term is by far the largest contributor to the friction factor 



for small pressure differentials, therefore making the approximations 



Equation(A-15) becomes 



i-- 



( |+ f x) 



where : 



(A-16) 



Kc = entrance coefficient 
Ke = exit coefficient 

whose values are determined from appropriate curves in Chapter IV, refer- 
ence (8). 

This is the equation used to determine the friction factor as 
presented in this thesis. 

If a further assumption is made such that 

■f * C|°« ££ £ <A-17) 



and equations (A-4) and (A-9) are substituted, Equation (A-17) becomes; 



~ P** A ^ \ f> 3 

V L ) <0 



(A-18) 
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Therefore, for a given matrix the friction factor is directly pro- 
portional to the cube of the porosity and indirectly proportional to the 
compactness . 

f - P 3 /jS (A- 19) 



COLBURN 1 FACTOR : 

We have A and maximum slope with which we get N £u from Table VIII 
or Figure 18 taken from reference (5). 

The j factor is defined as: 



j = N__ N 2/3 
J ST PR 



G c 



(A- 20) 



substituting (A-9) for G and multiplying by /A 



1 = ^ 3 — n 

^ mr A P 



2/3 



(A- 21) 



by definition: 



hA 

mCp 



Ar 2/T 

Therefore: j = N fcu — N pR (A-22) 

Note, from (A-l) Ac/A = r H/L, substituting f 3 / p for 

r R from (A-4) gives: 



2/3 



"|8 



(A-23) 



which leads to the relation that for a given matrix, j is directly pro- 
portional to the porosity and inversely proportional to the compactness. 

i ©c l °/|S 

HEAT TRANSFER POWER AND FLOW FRICTION POWER : 

The heat transfer power per unit area per degree temperature differ- 
ence is derived as follows: 
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j = N N, 

J ST 1 



2/3 



Solving for h: 



j = 



c p^ 



k 

c pM 



"r 

2/3 



p M 



2/3 



Now if c^, yU and N pR are evaluated at standard properties which by 
convention are taken as those for dry air at 500 °F and one atmosphere 

c = 0.2477 BTU/lb °F 
P 

/X = 0.0678 lb/hr ft 



1° 



= 0.0413 lb/ft 



Thus, the heat transfer power per unit area per degree temperature 



h STD = -° 2195 ( ”4^” ) (N r J > BTU/hr ft2 ° F 

The flow friction power per unit area is derived as follows: 



« ~ **■ & A? ( r -t) 



aP (t) = 



£l. l 
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£ -aP( £)v(£o) x f 






where £ x p V 

E = f 



4 (ek) S (f.*) ] 



Now if jj. and jO are evaluated at standard properties (previously 

defined) the flow friction power per unit area evaluated at the above 

standard conditions, becomes: 

’ STD’ 

e std = 1.11 X 10- 7 ( J -) 3 (\^) 3 f [hp/ft 2 J 

Figure 16 provides a comparison of the 20 ° skew and the "parallel 
plate" matrices from a heat transfer power vs. friction power point of 
view (hgjp vs. Eg^) . For comparison purposes the surface geometries 
have been adjusted to a common hydraulic diameter of 4r^ = 2 x 10 ^ft. 
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APPENDIX B. 



DESCRIPTION OF EQUIPMENT 

To obtain heat transfer data, a transient technique was employed. 
Briefly, this method consists of heating the test matrix, to a uniform 
temperature and then subjecting the matrix to a step change in fluid 
flow temperature. The fluid temperature downstream of the matrix is 
recorded versus time. 

Friction factor was calculated from pressure drop data obtained 
from static pressure taps located in the test section immediately up- 
stream and downstream of the matrix. 

The equipment necessary to complete the experiment falls into one of 
the following systems. 

1) Fluid source 

2) Flow metering system 

3) Temperature measuring system 

4) Pressure measuring system 

5) Fluid heater 

6) Matrix holder & test section casing 

FLUID SOURCE : 

The working fluid was, of course, air. Air was provided to the 
apparatus by placing the entire rig on the inlet to a lOOcfm, 6 stage 
centrifugal air compressor. The compressor requires an operating voltage 
of 12o/ at 90amp. 

FLOW METERING SYSTEM : (See Figure 1) 

The flow metering system consists of an orifice in a 3.08" inside 
diameter aluminum tube. The orifice section was constructed in accord- 
ance with the ASME specifications of reference (7) for standard D and 
D/2 pressure taps. A large flow range was obtained by means of 7 orifice 
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plates of throat diameter 2.310, 1.971, 1.232, 1.540, .462 and .308 
inches, respectively. Accurate control of the pressure drop across the 
orifice was maintained by a large sliding valve downstream of the ap- 
paratus on the inlet to the compressor and another sliding valve on the 
discharge side of the compressor. There is an additional valve on the 
inlet side of the compressor that is fully open or fully closed. It is 
normally opened at flow rates less than 150 lb/hr to prevent stalling the 
compressor. 

TEMPERATURE MEASURING SYSTEM : (See Figure 1) 

Without exception, all temperatures were measured with iron-constan- 
tan thermocouples. A single thermocouple measured the orifice tempera- 
ture (T q ) just upstream of the orifice at the inlet to the apparatus. 

This temperature was read in millivolts on a RUBICON CO. portable pre- 
cision potentiometer, converted to temperature and recorded for each run. 
Measuring is a group of 5 thermocouples in series, bound together and 
insulated from each other by teflon tape and inserted in a 1/8" diameter 
aluminum tube. The tube was placed in the air stream with a small window 
cut in it facing upstream. At this point ambient air temperature was 
measured. T 2 and T^ are 5 wire thermocouple grids permanently mounted in 
the test section cabinet. Their wires are connected in series to magnify 
the thermocouple output voltage. T^ is a movable 5 wire thermocouple 
grid located in the matrix holder and is placed adjacent to and down- 
stream of the matrix, it’s wires are also in series. It’s position is 
apparent in Figure 21. The thermocouples arranged in this manner provided 
flexibility and complete coverage of desired points of interest. If the 
thermocouple output voltages from T 2 and T^ are bucked against each other, 
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the temperature measuring system response curves may be generated on a 
Leeds & Northrup "Speedomax" temperature recorder. T2 vs gives the 
indication of temperature uniformity across the matrix and vs 
records the rate of change of the difference between upstream and down- 
stream temperatures with time. It is this record that provides the pri- 
mary data for the evaluation of heat transfer data. The precise relation 
of this record to the Colburn j factor is explained in Appendix A. The 
record is made on a Leeds and Northrup "Speedomax" Temperature Recorder. 

The recorder has variable chart speeds, which permits maintenance of accept- 
able cooling curves with regard to slope determination. The Speedomax also 
has a range of 0-20 mv. For this experiment, the recorder is precalibrated 
on a 3 mv. scale. 

There is one addition thermocouple grid that is not used in the experi- 
ment, but was used in checking out the apparatus. This grid would take the 
place of T^ and its purpose was to determine the temperature distribution 
across the cross section in both the vertical and horizontal directions. 

It was placed in a special holder that also held a velocity probe (See 
Figure 22) and was lined with balsa wood to provide a smooth and continu- 
ous internal flow passage. 

PRESSURE MEASURING SYSTEM : (See Figure 1) 

Pressure taps are located upstream and downstream of the orifice and 
upstream and downstream of the matrix. Each tap is connected by Saran 
tubing to an appropriate manometer or draft gage. The following instru- 
ments were used: 

1) Type "C" micromanometer, E. Vernon Hill & Co. 0-1,25" 

2) Ellison Differential Direct Draft Gage, 0-6", 0-8" 

3) Ellison Inclined Draft Gage, 0-6" 

4) Meriam Instrument Co. 0-60". 
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Any one or combination of these instruments can be used to measure 
the differential pressures of the orifice and the matrix. Frequent cross 
checks were made to ensure reliable operation. 

FLUID HEATER SYSTEM: (See Figure 1) 

The air heater section consists of 28 nichrome .0031" diameter wire 
heaters. The heater system is designed to heat 1000 pounds of air per 
hour, 20°F above ambient temperature. Nichrome was chosen because of 
its high resistivity, low thermal conductivity and low specific heat. 

The nichrome wire is wound on a bakelite frame, two heaters to a 
frame. The wires are 1/32" apart and there are 50 and 52 wires respective- 
ly, to each heater on a frame. The heaters on each frame are connected in 
parallel and then connected through a switch, to a variable voltage bus. 
There are two frames to a plastic holder, as shown in Figures 23 and 24. 

In these figures, the heater frame and its electrical connections can be 
clearly seen. 

A schematic wiring diagram for one frame is shown below. All 14 
frames are wired in this manner, putting all 28 heaters in parallel with 
the capability of switching on or off two at a time. 
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The 28 heaters arranged in this way, provide a variety of flexi- 
bility. The number of heaters in use, should decrease as the flow rate 
decreases. Voltage variation can be obtained by a GENERAL RADIO COMPANY 
"VARIAC", 0-240 V , 8a, 50-60 cycles. It is recommended that a voltage of 
220 V or in this vicinity, be maintained, at all flow rates, and the 
number of heaters used, decreased until the number of heaters is four. 

At this point, the voltage should be decreased, as necessary, to obtain 
a 20° temperature rise, the number of heaters remaining constant. A 
note of caution: Do not energize a heater, unless the air is flowing, 

otherwise, the wires will sag due to thermal expansion and short circuit- 
ing will occur. 

MATRIX HOLDER AND TEST SECTION : 

The matrix holder and test section casing are made from polyethylene 
plastic. This material was chosen because of its low heat capacity. 

The test section casing holds thermocouples and T^ permanently 
mounted upstream and downstream respectively. It also contains the up- 
stream and downstream static pressure taps used in determining the fric- 
tion factor. The parts making up the casing were machined to close toler- 
ances, to ensure a tight fitting matrix holder, and a minimum of air leaks. 

The parts comprising the matrix holder, as implied above, were machin- 
ed to close tolerances to ensure a tight fit. Great care was taken to 
ensure smooth alighment of the flow channel through the heaters and the 
test section. The size of the matrix holder was arbitrarily chosen. 

The flow channel dictated the width, and the length is such that matrices 
of various flow lengths, up to 3", may be tested. The flow channel is 
3-3/16” x 3-3/16”. 

The matrix holder has one special provision and that is for the movable 
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thermocouple T^* A small hole is drilled in the front, and the wires 
are led through. The matrices are placed in the holder and held in 
place by bits of Styrofoam plastic cut to the size necessary to hold 
the matrix square in the flow channel. 

That completes a description of the various systems involved in 
the apparatus, but there is one other item that bears mentioning. Just 
upstream of the heater section, is a section provided for screens to 
straighten the velocity profile. Provisions were made for 1 to A screens. 
One 60 mesh stainless steel screen was found to be sufficient, to straight- 
en the flow at the flow rates used (40-1000 lb/hr.) 
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APPENDIX C 



CALCULATION OF GEOMETRIC CONSTANTS 
FOR PERFORATED MATERIAL 

As demonstrated in Appendix A, the surface geometry has a strong 
influence on the heat transfer and friction factor curves. Specifically, 
the effects of p and jj were examined. 

In using perforated materials, the presence of the slots, require 
modification of these two parameters. It is the intent of this Appendix, 
to demonstrate the calculation procedure utilized to determine the geo- 
metric parameters, on which the graphical and tabular results for the 
matrices constructed from perforated material are based. 

The material used in the perforated plate matrices, was a pure nickel 
type 160/40TV, slotted opening plate material, .0022 inches thick. This 
material is a product of Perforated Products, Inc. 

The basic geometry of the plate is: 



\ JWS, 



^ 622 cm 

‘Y~ 



~JMK "/////A WE. 



"T 



/////A 



Using the average between the top and bottom of the tapered slots, 
consider an average slot pattern as: 
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CROSS HATCHED AREA A , = (.063) (.016) = .001008 cm 2 
cn 

AREA OF SLOT IN A^ = (. 040) (. 0020) = .00008 cm 2 

Now define Y = Aslot = AREA REDUCTION DUE TO SLOTS 



Y- 



.00008 

.001008 



7.937. 



Therefore, plate Porosity (p g ) = (1 - V ) = .9207 



Now define 



£ 



CROSS SECTION AREA OF SLOTS 
CROSS SECTION AREA 



CONDUCTION AREA 
REDUCTION DUE TO 
SLOTS 



6 = 



(.040) (.0022 ) 
(.063) (.0022) 



= 63.57. 



The following calculations are for the 20° skew perforated plate 



matrix: 



DIMENSIONS OF PLATE: 


4.285" x 1.98" x .0022' 


NUMBER OF PLATES: 


142 


DIMENSIONS OF MATRIX: 


3.255" x 3.248" x 1.98' 


WEIGHT OF MATRIX: 


275.5735 grams 


MATERIAL CONSTANTS; 


fl = 555 lb / ft 3 




V 38 ‘ 7 BTU/ hr°F ft 




V * 1065 BTU/ lb °F 



PLANE SURFACE AREA = H A D (*'-*%*) V • 98 ) ( 2 ) = 16.733 ft 2 

144 

TOTAL HEAT TRANSFER AREA (A) = PLANE SURFACE AREA - AREA OF SLOTS 
= (1- ^ ) SURFACE AREA 
= (.9207) (16. 733) 



FRONTAL AREA (A fR ) 



A = 15.41 ft 



3.255 x 3.248 = 10.572 



. 2 
m 
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FREE FLOW AREA = FRONTAL AREA - CROSS SECTIONAL AREA CF METAL 

2 2 
= 10.572 in - 4.285 x .0022 x 142 = 9.233 in 



A = .06416 ft 2 

c 

CONDUCTION AREA (A g ) = (142) (4. 285) (.0022) (1 - £ ) 

= (142)(4.285)(.0022)(.375) = .50198 in 2 

A = .003486 ft 2 

s 



MAXTIX VOLUME (V ) = 3.255" x 3.248" x 1.98" = 20.933 in 3 

MATERIAL VOLUME (V g ) = 142x4. 285"xl . 98"x. 0022" = 2.650 in 3 

WEIGHT OF MATERIAL (W g ) = 275.5735 gms./453.6 gm/lb = .6075 lb. 

MATERIAL VOLUME (\T) - W / n = .6075 lb/555 lb/ft 3 = .001095 ft 3 
s s Ps 

The porosity of a matrix has been defined as: 



MATRIX FLOW VOID VOLUME 

MATRIX VOLUME m ~ s 

Vm 



Vs 

V 

m 



P 



. 001534 

.01211 



.8733 



P 



V s . . 001095 _ 

V = .01211 " 1 



.09042 = .90958 



The compactness of a matrix has been defined as 



P 



TOTAL HEAT TRAN S FE R AREA 
MATRIX VOLUME 



Here, 



we will redefine 






to include the effect of area reduction: 



P- 



(PLANE SURFACE AREA) (1 - V ) 
MATRIX VOLUME 
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For Y = 0 



P- 



(16-733H1) 

.01211 



For 



/ = .0793 

P- 



(16,733) (.9207 ) 

.01211 



1272.18 ft 1 ft 3 



The hydraulic diameter is conventionally defined as: 



4 FLOW CROSS SECTION AREA 
°H = WETTED PERIMETER 

Appendix A shows that hydraulic diameter can also be expressed as: 



d h ■ 4 p/ ( s 

To be consistent with the conventional definition, for the case of 
a perforated material, is used vice , since the slots have no 

effect on the hydraulic diameter. 

Finally then, 



. 8733 

1381.75 



.002528 ft 
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APPENDIX D 



SAMPLE CALCULATIONS 



The following sample calculations are based on recorded data from 



Run 43 for 20° skew matrix constructed from perforated nickel sheets. 
This is chosen primarily because the basic geometric parameters are com- 



puted in Appendix C. 



RECORDED DATA : 



D 


= 


2.310" 




o 








> 

<-a 

o 


= 


2,61 in. 


h 2 o 


AP 


= 


3.96 in. 


H o 0 


m 






2 


P 


_ 


4.69 in. 


H_0 


s 






2 


T 


= 


68.0 °F 





MATERIAL CONSTANTS : 

c g = .1065 BTU/lb °F 

k g = 38.7 BTU/hr ft°F 

PRECALCULATED GEOMETRIC PARAMETERS : 

A = 15.41 ft 2 

A. = .0734 ft 2 
fr 

A = .06416 ft 2 
c 

X = .003486 ft 2 
s 



P „ = 29.963 in. Hg 

atm 

CHART SPEED =4.0 in/ sec 
W s. 6075 lb 



(Re: Appendix C) 

jO = .8733 

|3 = 1272.18 ft 2 /ft 3 
r R = .000632 ft 



L = 1.98 inches ( — denotes corrected for slot effects) 
DETERMINATION OF MASS FLOW : as outlined in Reference (7) 

m = 359 K D 2 F Y J A P T” 
o a v o 

where: Y * 1.0 Fig. 40b Ref. (7) 

F g « 1.0 Fig. 38 Ref. (7) 
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m 



(359) (2. 31) 2 K ''I (2.61) ( (.07522) 



m = 849.40 K 



assume K = .749 




Then m = (849.40) (.749) = 636.21 lb/hr 

Reynold's number for the pipe N = 

Rp 





= 7.35 x 10 



From Table 5, Ref. (7) /9= 0.75, N R = 7.35 x 10 4 , K = .749 



Therefore 



m = 636.21 lb /hr 



MATRIX REYNOLD'S NUMBERS 

There are actually two Reynold's numbers to be calculated for the 
matrix. One is for isothermal flow and is associated with the .friction 
factor, the other is associated with heat transfer. They differ only in 
the temperature at which the absolute viscosity JU. is evaluated. ^ 

(friction) is evaluated at ambient temperature (T q ) . JA ^ (heat transfer) 
is evaluated at T q + 10°F (an estimate of the average bulk fluid tempera- 
ture) . 




(^(^ff^) ft(636.21 lb hr 



= 618.58 



(15. 41ft*) (.04405 lb/ft hr) 




CONDUCTION PARAMETER 




k A 



.00536 
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MAXIMUM SLOPE 



From Figure 26 the slope of the cooling curve was determined as 
described in Appendix A. The slopes of three curves were determined and 
the average of the three was used in the calculations. 



MAXIMUM SLOPE = 



C s / dT_\ _ C s 



SLOPE xCHART SPEED 



d©/max 



(.6075 lb) (. 1065 BiU /lb°F) = .06470 B1U /°F 



C = m c p = (636. 21 lb /hr) (. 24 BTU /lb°F) 



.65745 



SLOPE = ,16014 in (See Figure 26) 

Finally: 

MAXIMUM SLOPE = > <‘ 10775 > 4 > 

COLBURN j FACTOR 

With X and MAXIMUM SLOPE enter Table VIII and find N = 4.40. 

tu 

Linear extrapolation will give desired accuracy. 



Now, 



n pr 



2/3 



from Appendix A, Eq. (A-22) 



therefore, j = * 15^41 (.796) (4.40) = 

FANNING FRICTION FACTOR 



.01458 



As derived in Appendix A Equation (A-16) the friction factor is: 



f . 







L 



AP = 3.96 in H„0 x 5.204 
m 1 



psf 

in H 2 0 



= 20.61 psf = .1431 psi 
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Af = 4.69 in H o 0 x 5.204 P 8 - - -- = 24.41 psf = .1695 psi 

s 4 in 

P = 29.963 in H x .4912 P sl „ = 14.7178 psi 
atm g in 

P = P - P = 14.7178 - .1695 = 14.5483 psi 
1 atm s 



P = = 14.4768 psi 

mean 2 

_ ( U„> (144) = (14.477) (144) = lb 3 

r mean RT (53.3) (528) ' U/4U/ ' 

I m 

m 636.21 lb/hr « ,, „ ,..2 „ , lb , ^ 2 

G = ~T~ = . 06416 " ft 2 ' = 9915 * 9 lb/hr ft = 2 * 754 /sec ft 



.000632 
r H/L = 1.98/12 



= .00383 



K = -.738 



K = .949 from Chapter IV Reference (8) 



(1 + p) = 1 + (.8733) = 1.7626 

AP 






P 



mean 



_ ( 64.4) (.07407) (20.61 )= 12.9579 

,2 " 7.587 



AP 2. 

_E C 1+ P ) 



P Sl . (1.761) = .01742 

14.4768 psi v ' 



f = (12.9579 - .01740 - . 211) (. 00383) = (12 . 7295) (. 00383) 

= .04859 

All thermal properties ( , Npp,k g , and C g ) jere obtained from 

References (12) and (13). The temperature at which the properties were 
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evaluated was the ambient temperature with the exception of the absolute 
viscosity for heat transfer ( )• 
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APPENDIX E 



A FORTRAN COMPUTER PROGRAM 
FOR DATA REDUCTION 

A short computer program was devised. Figure 27, so that from basic 
experimental data, the end result is: 

1) Conduction parameter (A) 

2) Maximum slope 

3) Reynold's number for heat transfer 

4) Reynold's number for friction 

5) Fanning friction factor 

It was anticipated that this program along with a subroutine would 
ultimately produce N fc ^ and subsequently the Colburn j factor. 

Time did not permit the adaption of Moreland's work (11) to a sub- 
routine but the program is still useful. It's value when data for a 
single matrix is to be reduced is nebulous. However, the reduction of 
data for multiple matrices can be accomplished with a considerable savings 
of time. 

The following is a list of inputs necessary in Fortran and English 
language along with the required dimensions: 

CARD NO. 



XML 


Matrix length 


ft 


AS 


Solid conduction area 


ft 2 


RH 


Hydraulic radius 


ft 


POR 


Porosity 


- 


AFL 


Flow cross section area 


ft 2 


AHT 


Matrix heat transfer area 


ft 2 


WM 


Matrix weight 


lbs 
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CARD NC. 



40 EK Exit flow coefficient 

CK Entrance flow coefficient 

SK Solid thermal conductivity 

CM Matrix material specific heat 

N Number of runs 

10 DO Diameter of orifice plate opening 

DELPO Pressure drop across orifice 

GAM Specific weight of fluid 

ATMP Atmospheric pressure 

COR Orifice correction K 

UH Fluid viscosity for heat transfer 

UFR Fluid viscosity for friction 

10+2 CS Chart speed 

DELH Pressure drop across matrix 

HS Static pressure matrix inlet 

SLO Slope of generalized cooling curve 

TEMPO Orifice temperature 

The output of this program in Fortran and English 
their dimensions is: 

EMDOT Mass flow rate (m) 

RNUMP Reynold’s number for the pipe 

RNUMH Reynold's number for heat transfer 

CAPFL Thermal capacitance of the fluid (C) 

CAP: Thermal capacitance of the solid (C g ) 

CONDPAR Conduction parameter ( X ) 



BTU/hr °F ft 
BTU/lb°F 



inches 
inches I^ 0 
lbs/ ft 3 
inches Hg 

lb/hr ft 
lb/hr ft 
sec/ inch 
inches I^O 
inches HjO 
(inches) * 

°F 

language along with 
lb/hr 



BTU/sec °F 
BTU/ °F 
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SLOMAX 



Maximum slope of generalized cooling curve 



G 



Mass flow velocity 



lb/hr ft' 



2 



RNUMFR 



Reynold's number for friction 



FFR 



Fanning Friction factor 



Reduction of data is not complete since j remains to be found. 

With SLOMAX and CONDPAR, enter Table VIII, and find N . The Colburn j 
* ' tu J 

factor follows from equation (A-22), Appendix A. 

The program is flexible in that an unlimited number of matrices 
may be reduced by changing the first card of the program appropriately. 
(L = the number of matrices to be reduced.) 
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